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Abstract 

The method of smoothed particle hydrodynamics (SPH) is apphed for ultra- 
relativistic heavy-ion collisions. The SPH method has several advantages in study- 
ing event-by-event fluctuations, which attract much attention in looking for quark 
gluon plasma (QGP) formation, because it gives a rather simple scheme for solving 
hydrodynamical equations. Using initial conditions for Au+Au collisions at RHIC 
energy produced by NeXus event generator, we solve the hydrodynamical equation 
in event-by-event basis and study the fluctuations of hadronic observables such as 
dN/dy due to the initial conditions. In particular, fluctuations of elliptic flow coef- 
ficient V2 is investigated for both the cases, with and without QGP formation. This 
can be used as an additional test of QGP formation. 



1 Introduction 

One of the central issues of the high-energy heavy-ion physics is to investigate the phases of 
the hadronic and quark matter. Hydrodynamical descriptions|l|], which have a rather long 
history, may be a powerful tool for studying such phase diagrams because we can handle 
the equation of state directly in the model. However, to extract precisely information 
about the equation of state, we must be careful in comparing the model predictions 
with experimental data because procedures of event averaging may make signal of QGP 
ambiguous. 

Event-by-event analysis [§] is one of promising ways to extract clear information from 
the experimental data, in particular at RHIC and LHC energy regions. The fluctua- 
tions at the initial stage (for example, the initial energy density fluctuation) must affect 
the formation of QGP and the later space-time evolution of the whole system. Fig.l 
shows the energy density (counter plot, at z = fm plane) for a Au+Au collision at 
y/s = 200 A GeV, 6 = fm, produced by NeXus code|^. As shown in Fig.l, the fluctua- 
tions at the initial stage of relativistic heavy-ion collisions are not negligible. To achieve 
an event-by-event analysis using relativistic hydrodynamical model, the numerical code 
must deal with arbitrary initial conditions and equation of state with suitable calculational 
speed and precision. The method of smoothed particle hydrodynamics (SPH)^, which 



was first introduced for astropliysical applications 0, satisfy such requirements. The main 
characteristic of SPH is the introduction of "particles" attached to some conserved quan- 
tity which are described in terms of discrete Lagrangian coordinates. This feature of the 
model enables to easily carry out calculations of the ultra-relativistic heavy-ion collisions, 
which is accompanied by prominent longitudinal expansion. Another advantage of using 
SPH method is that we can choose a suitable precision in solving hydrodynamical equa- 
tions because of the smoothing kernel, which allows to smooth out often unnecessary very 
precise local aspects and then achieve a reduction of computational time. This is a very 
profitable point in the event-by-event analysis. 

In the following section, we briefly formulate the entropy-based SPH, using the vari- 
ational approach]^. In Sec. 3, to check the performance of SPH, we apply it for solving 
problems of relativistic hydrodynamics whose solutions are known analytically or numer- 
ically. Then, in Sec. 4, we demonstrate hydrodynamical evolutions of both resonance gas 
and quark gluon plasma starting from the same simple initial conditions. In Sec. 5, af- 
ter giving a brief explanation how to take the initial conditions of NeXus into our code, 
we present typical hydrodynamical evolution of an event at RHIC energy. To introduce 
freeze-out process into our code, we formulate Cooper-Frye formula in terms of SPH and 
show several results on single-particle spectrum in Sec. 6. Then, we study the event-by- 
event fluctuation of elliptic flow coefficient f20 at RHIC energy P] in Sec. 7, showing that 
the dispersion is large and it does depend on the equation of state used. Finally we close 
our discussions with some concluding remarks. 



2 Relativistic SPH equations 
2.1 SPH representation 

To specify the hydrodynamical state of a fluid, we need to locally know the thermody- 
namical quantities such as the energy density e and the collective velocity fleld, besides 
equation of state of the fluid. If local thermal equilibrium is well established, the ther- 
modynamical variables may be expressed by smooth functions of the coordinates. Let us 
choose space points to specify the state of the whole system. To know the physical 
quantities in other space points, we may use an interpolation using the known N points. 
(Hereafter we use the word 'particles' instead of 'points'.) 

In SPH, for an arbitrary extensive thermodynamical quantity of A, the density (in the 
space-flxed frame) a* is parametrized in the following way: 

N 

a*{x,t) = Y,u, W{x-x,{ty,h) , (1) 

i 

where W{x — Xi{t); h) is a positive deflnite kernel function with the properties 

J dx' W{x' — x) = 1 and W{x' — x,h) = 5{x' — x) when h ^ , 
where /i is a smoothing-scale parameter. Using the normalization above, we have 

N 

Atotai = J d^x a*{x, t) = Y,Ui . (2) 
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From eq.(H) the quantity Vi can be interpreted as a portion of A which is carried by the 
ith. particle. Its velocity, identified as the velocity of the fiuid at Xi{t)^ satisfies 

vm = ^ . (3) 

If A is a conserved quantity, z/j must be constant in time. In this case, from eqs.(|lD and 
d), we obtain 

^a*{x,t) = -V -jix^t) , (4) 

where 

3{X, t)=Y. VV{X - X,{ty, h) (5) 

i 

is interpreted as the current density of the thermodynamical quantity A and eq. @) is then 
the corresponding continuity equation. 



2.2 Equations of motion for SPH particles 

To know the time evolution of the fiuid, we have to find the equation of motion for 
each SPH particle. It can be obtained by the variational method @] with the following 
Lagrangian[Q 

Ls^^{{xi, V,}) = -Y.^i ^ (6) 

i 

It should be noted here that the factor z/j/a* can be regarded as the volume Vi (in the 
space-fixed frame) occupied by the ith 'particle'. Then each term in eq.(|^) is equal to 
Ei/'^i, where Ei and 7^ are the energy (in the rest frame) and the gamma factor of the 
ith particle, respectively. Then we also write 

The equations of motion of SPH particles are obtained minimizing the action JgpH with 
respect to {xi} with the constraint given by eq.(^) taken into account 



6Ispu{{xi, Vi}) = -5 (dt ^ = . (8) 
J li 

Using the adiabatic condition 

= P6V, = 6{^) , (9) 

we obtain 

d P. -^ P. r P. R 1 

VW{xi - Xj] h)=0 , (10) 



d . Pi + Si 

tS^^ — 

dt a,- 



— [Vi '-liVi) + 

3 



P P 

*2 ' *2 



where Pi is the pressure. Because of the adiabatic condition, eqs.(|T3) are the equations 
of motion for the ideal fiuid. For the extension to the non-adiabatic case, see ref.iffl. 
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2.3 Entropy based SPH in general coordinate system 

As a possible candidate for the thermodynamical extensive quantity a*, we may consider 
baryon number density, entropy density, energy density and so on. In the case of energy 
density, we have (cf. eq.(|D and (|^)) 

z/, = E,{x,,v,) . (11) 

However, this is only the internal energy, so I'i is not constant with respect to the variation 
of Xi in eq.(H), introducing an additional complication without any practical merit. On the 
other hand, for entropy or baryon number density, for example, z/j can be kept constant. 
Since we will apply SPH to ultra-relativistic heavy-ion collisions, baryon number density 
is not so suitable choice because it is expected to be very small in the central rapidity 
region. Let us choose the entropy density for a* in eq.(|l|). 

In relativistic heavy-ion collisions, initial conditions of fluid are often given at constant 
r = yW—z^. Then, it is convenient to formulate SPH in the coordinates 

= T = \Jt'^ — 2;^, x"^ = X, = y, = ri = - In . (12) 

2 t — z 

The SPH equations in generalized coordinate system are easily derived in a similar way|^. 
In the case of the coordinate above, the SPH equations of motion should be written as 

s(r, a;) = ^ V z/iiy (ic - ajj : /i) , (13) 



Pi Pi 



9 o2 9 o2 



VW[xi-xf,h] = {] . (14) 



3 Numerical check of SPH 

To check the performance of SPH formulated in Sec. 2, we applied it to a problem of 
relativistic hydrodynamics whose solution is known analytically. When the longitudinal 
rapidity a and the transverse one {3 are given as 

1 r + y/g;^ + 

a = V, /^=oln , . ^ , (15) 

/ r — va^ + y 



it is known that the entropy density is given by 

So 



(16) 



where Sq is a constant which is determined by the initial conditions. A comparison of the 
numerical results by SPH with the analytical solution eq.(|16D is shown in Fig. 2, where 
the equation of state used was = 1/3. It is verified that SPH solves this problem 
correctly. For a further check, we applied it also to another problem, namely, relativistic 
hydrodynamics with the longitudinal scaling and the Landau-type initial condition in the 
transverse directions, 

a = 7], p = 0. (17) 

Figure 3 shows a comparison of the SPH results with those obtained by Hama and 
Pottag[§]. The initial conditions and the equation of state used are exactly the same 
as in ref. [P[. We found that SPH works well also in this problem. 
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4 Equations of state and hydro dynamical evolutions 



To complete the theory of relativistic hydrodynamics, the equation of state (EoS) for the 
fluid is required. We use two possible EoS in the present study, namely, a resonance 
gas EoS|jlO| and an EoS containing both resonance gas and QGP via a first order phase 



transition [Pm. The EoS are parametrized as follows (See also Fig.4.): 

P/s = 0.20 ( constant for whole e region ), (18) 

for the resonance gas (RG) and 

r 0.20 0.00 GeV/fm^ < e < 0.28 GeV/fm^ = ei 

P/e = I 0.056/e 0.28 GeV/fm^ < e < 1.45 GeV/fm^ = €2 (19) 

[ l/3-AB/3e 1.45 Gel// /m^ < e 

for the quark gluon plasma (QGP) EoS. The bag constant B = 0.32GeV/ fm^ is deter- 
mined to give the critical temperature = 160 MeV. A remark should be made here that 
when we call QGP EoS, here and in the following sections, we are actually considering all 
these three states according to the value of e. 

We show some results of the hydrodynamical evolution using such RG and QGP EoS in 
Figs. 5 and 6. For the both cases, the initial (tq = 1.0 fm) entropy density s, longitudinal 
and transverse rapidity a, (3 are set as 

s = 40[l//m^], a = 7], /3 = . for ^/x' + < 7.0 /m, |r/| < 7.0 . (20) 

As seen in Figs. 5 and 6, the space-time evolution of the fluid strongly depends on the 
EoS. In comparison with the RG case, the temperature of the fluid in QGP EoS drops 
much slowly due to the existence of the mixed phase, where the pressure is constant. 



5 Initial conditions produced by NeXus 

The most fundamental assumption of hydrodynamical models is the thermalization of 
partonic or hadronic systems produced after high-energy heavy-ion collisions. In our 
study, we set the initial conditions using information from Nexus event generator [Q, in 
the following 2 steps. 

1. Align all primary hadrons on the r = tq hypersurface. Here we assume that the 
primary hadrons move freely with the momentum p^. In the present study, we use 
To = 1.0 fm. 

2. The initial energy density and momentum density are estimated using the interpo- 
lation kernel W . 

e{TQ] x) = J2p'^iW{To, X - Xi] h) , 

i 

p{to; x) =J2 Pi^ i^o, X - Xi] h) , 

i 

where Xi = {xi^yi^rji) are the coordinates of the i-th primary hadron after the 
alignment on r = tq. The smoothing scale parameter h = 1 fm is used, which 
is the rough size of the hadrons. The collective velocity field is then given by 
v{to,x) =p{To,x)/e{To,x) . 
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In Fig. 7, we show the hydrodynamical evolution of a typical event of Au + Au collision 
at energy ^/s = 130A GeV for the impact parameter b = 7.0/m given by our hydro code, 
SPheRIOQ, connected to the NeXus event generator. 



6 Cooper-Frye formula in SPheRIO 



The Cooper-Frye formula|T2| is widely used in calculations of single particle spectra using 
hydrodynamical models because of its simplicity. So, as the first trial to evaluating single 
particle densities, we shall follow it. 

In SPH, using eq.([TB|), the usual Cooper-Frye formula can be rewritten as 



E- 
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E 



F.O 



(T 1 

TdxdydTj 3^ W{x - x^ir); h) ^^^^^./y^ ^ 



where is defined by 



0",, 



Tlx ^rj 
Tt/^ Tti^ T'lj'T' 



[21) 



(22) 



and n,, 



{nT-,nx,ny,nrj) is the normal vector of the isothermals: 



n. 



dT dT dT dT 

dr' dx^ dy^ dt] 



(23) 



The integration is usually done over a constant temperature hypersurface, T = Tj. When 
the spline-type kernel functionj^, |^ is used, W{x — x'; h) = for \x — x'\ > 2h. Then, 
within the approximation of small h, 



E' 



dp^ 



9 ^ 



(27r) 



E 



/ rdxdydrj Wix — XjirJx)) ; h) 

j J T su^ 

afj,W{x — Xj ; h) 



dxdydrj 



sv>- 



(24) 



J Ut 



where Xj gives the point where j-th SPH-particle crosses the freeze-out hypersurface. 
The result of the single-particle spectra, l/nirj^dN / dm^, dN/dy and dN/dip for Au+Au 
collisions at y/s = 130A GeV, b = 7.0 fm are shown in Figs. 8, 9 and 10, respectively. 
The freeze-out temperature used is Tf = 140 MeV. 



7 Event-by-event analysis of hadronic observables 
7.1 Elliptic flow coefficient V2 



The flow phenomena ||13|| can be important candidates of QGP signals because they may 
carry much information about EoS during the expansion of quark or hadronic matter. 
In particular, the elliptic flow coefficients is one of the interesting observables since it 
is sensitive to the early stage of high-energy heavy-ion collisions. Because it is expected 
that different EoS responds to the initial fluctuations in different way, the measure of 
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fluctuation, for example, Sv2 = y (^1) ~ ('^2)^ is also very important. For this purpose, 
we investigate V2 distribution using different EoS. In Figs. 11, 12 and Table 1, we show 
results of the V2 distribution for Au + Au collisions at energy ^/s — ISO^l GeV in RG and 
QGP EoS cases. 

As seen in Table 1, the fluctuations of flow coefficient V2 in QGP is about 20 ~ 30% 
smaller than that in RG (cf. the difference of (^2) in RG and QGP. It is about a few 
%.). The h dependence of 5v2 is also found. Hence, it is expected that 5v2 also brings us 
useful information to infer on the phase passed in the early stage of high-energy nuclear 
collisions. 

7.2 Fluctuation of the slope parameter 

It is also interesting to study the fluctuation of the so-called slope parameter in the trans- 
verse momentum spectra. The slope parameter T can be obtained by fitting 1 / rrir^dN / diriT- 
with a function A/"exp [— m-p/T]. In the present study, < — rrij^ < 1.0 GeV interval 
is used for the fitting. The results are shown in Table 2. We found T and 5T of QGP is 
systematically smaller than that of RG. 

7.3 Multiplicity fluctuation in the central region 

Other interesting measurement may be multiplicity fiuctuation. Let us define the fiuctu- 
ation measure as 



We investigate the EoS dependence of the multiplicity fiuctuation in two central {y — 
0) rapidity windows Ay and as function of the impact parameter h. The results are 
summarized in Table 3. As seen there, the multiplicities and their fluctuations (5n have 
a clear EoS dependence. See Figs. 13 and 14. When the initial energy density is bellow 
82 the multiplicity fluctuation with QGP EoS is expected to be larger than that of RG 
without flrst order phase transition. This is due to the existence of the mixed phase (see 
Appendix). The values of ratio bnj (n) of QGP is shghtly larger than those of RG. 

8 Concluding remarks and discussions 

Wc developed a new hydrodynamical code, SPheRIO, based on the smoothed particle 
hydrodynamics (SPH), for studying high-energy heavy-ion collisions in event-by-event 
basis. In this study, we set the initial conditions based on the results of NcXus event 
generator and used two types of possible equation of states (EoS), ie., the resonance 
gas (RG) and quark gluon plasma (QGP) EoS. In the freeze-out process, single particle 
spectra are calculated by the Cooper-Prye formula for each set of initial conditions and 
EoS. Analyzing these single particle spectra, the fluctuations in flow coefficients, slope 
parameter T and multiplicities in the central rapidity regions are investigated at RICH 
energy. 




(25) 



where 




(26) 
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As seen in the results of Table 1, we found remarkable difference in 5v2 between RG 
and QGP EoS predictions. This means that the measure of the fluctuations is also helpful 
to infer the phase passed in the high-energy heavy-ion collisions. In the QGP EoS case, 
the fluctuations are about 20 ~ 30 % smaller than those in RG case. Our predictions of 
the average value of V2 is shown in Fig. 15, compared with STAR Collaboration data0. 
It is seen that the main trend of the data is reproduced, but the fluctuation is quite large. 

The fluctuation of the slope parameter T is small, a few % of T. We found the 
fluctuation of multiplicity has a clear EoS dependence. The impact parameter dependence 
of 5n{y, Ay) can be a good measure for the detection of the first order phase transition. 

Acknowledgments: The authors acknowledge stimulating discussions with K. Werner, 
O. Socolwski Jr and H.J. Drescher. This work was partially supported by FAPESP (contract 
no.s 98/02249-4 and 98/00317-2), PRONEX(contract no. 41.96.9886.00), PAPER J (contract 
no. E-26/150.942/99) and CNPq. 



A Entropy fluctuations in the first order phase tran- 
sition 

A.l Mixed phase dominant case: 

If in the most part of the fluid the energy density at the initial time Tq varies only in the 
interval ei < e{x) < 62, where ei and 62 are the critical densities of the RG and QGP 
phases, respectively, the produced entropy S of the system is 



S 



d^x si 



e{x) + 5e{x) 



£1 



d^x Si 



e[x] 



1 



/5e{x) 
^Kx) 



(Al) 



where = O(constant) is the sound velocity in the mixed phase(Mix) and Si is the 
critical entropy density corresponding to ei. e{x) and 5e{x) denote the event averaged 
energy density and fluctuation, respectively. Here, we assume that 5e{x) <^ s{x), so that 
the mixed phase is dominant in all the events. Then, the dispersion of the total entropy 
S is 



(A2) 



In the case the phase transition is absent, so that the fluid is described by RG EoS, the 
same events with the same energy distribution and fluctuation gives 



d^x Si 
d^x Si(^ 



e{x) + 5e{x) 



e{x)\ol% 



5 / fe(x) 
6 V six) 



5 /fe(x) 



and 



72 V six) 
e(a;)\5/3 25 / /5e{x)\2' 



+ 



e[x] 



(A3) 



(A4) 
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respectively. We used the sound velocity Cr = l/-\/5( const ant) here. So, the appearance 
of the mixed phase produces 



Hence, 



(A5) 



A. 2 QGP phase dominant case: 

If the most part of the fluid is in the QGP phase at initial time Tq, the total produced 
entropy S of the system is 



S, 



QGP 



e{x) + 6e{x) 
e2\ /e(x)\3/i 



d'x .i(^) C-^Y'^ \ 1 + ^(^) - :^r^Y + 



3 i'6e{x)\ 3 /6e{x)-^2 

4 V e{x) ^ " 32 V £[x) 



and 



^QGP 



d X s 



2/£2\2 /g(x)x3/2 9 // fe(x) x2 
ApJ V J 1fi\ V pr-r^ / 



^Sl^ ^ 82 ' 16 \ ^ 

where = l/v^(constant) is assumed. Then, we have 



D 



QGP -'^RG 



3 /^^2\2 /e(x)\3/2 25 /£:(x)\^5/3 



^^^^^^ leii? 



(A6) 



(A7) 



be[x)\2 
e{x) 



From this relation, we can conclude that when the energy density e{x) becomes larger 
than a certain value, (0.81)^ {82l^\f'^\ ~ H GeV/ frn? for the most part of the fluid, 
the dispersion of the QGP case becomes smaller than RG case. 
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Figure 1: Initial energy density (counter plot z = fm plane) of a typical event of 
Au+Au collision at energy ^/s = 200 A GeV, impact parameter b = Ofm, produced by 
NeXus event generator 0. 
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Figure 2: Comparison of results by SPH (open circles) with the analytical solution, eq. ([l6D . 
The sizes of the parameters h-j:{= hy), h^^ and dr used were 0.030, 0.067 and 0.05, respec- 
tively. 




Figure 3: Numerical results by SPH (open circles) for the problem of longitudinally scaling 
hydrodynamics using the Landau-type initial conditions in the transverse directions. The 
solid line has been obtained numerically by Hama and Pottag [^. 
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Figure 4: Equation of state which contains the QGP phase transit ion |TT|. For resonance 
gas case, c„ =0.2(constant) is used. 
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RG, tau = 1 JOim 




Figure 5: Hydrodynamical evolution (temperature T [GeV] aX z = Ofm plane) using RG 
EoS. The smoothing scale parameter h used in x-y space is 1.0 fm and in rj space 1.0, 
respectively. 
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QGP, tau =2.0 lim 




Figure 6: Hydrodynamical evolution (temperature T [Gel^] at z = Ofm plane) using 
QGP EoS. At r = 2.0 fm, we can observe QGP phase (above T > Tc = 0.16 GeV) and 
mixed pliase(T = Tc) clearly. At around r ~ 3 fm, QGP phase almost disappear. Mixed 
phase survives up to about 25 fm. ^ 



Au+Au 130A GeV, b=7.0 fm, iau=1 JD ifm 




Figure 7: Hydrodynamical evolution (temperature T [GeV^]) of a typical event in Au+Au 
collision at energy ^/s = 130A GeV for the impact parameter b = 7.0 fm {at z = fm, 
impact parameter direction is parallel to x direction). The initial conditions are produced 
by NeXusg. The EoS used is QGP one. 
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Figure 8: mT-spectra in Au + Au collisions at ^/s — 130A GeV, b — 7.0 /m, for 5 events 
(open circles for RG and filled circles for QGP). 
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Figure 9: dN/dy in Au + Au collisions at a/s — 130 A GeV, b — 7.0 /m, for 5 events (open 
circles for RG and filled circles for QGP). 
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Figure 10: dN/d(p in Au + Au collisions at ^/s = 130A GeV, b — 7.0 fm, for 5 events 
(open circles for RG and filled circles for QGP). 
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Figure 11: The distribution of elliptic-flow coefficients V2 for Au + Au collisions at energy 
y/^ = ISOA GeV{b = 10 fm). RG EoS is used. 
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Figure 12: The distribution of eUiptic-flow coefficients V2 for Au + Au collisions at energy 
= 130A GeV{b = 10 fm). QGP EoS is used. 
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Figure 13: EoS dependence of {n{y,Ay)) {y — 0, Ay = 1.875) as function of impact 
parameter b. 




Figure 14: EoS dependence of 5n{y,Ay) [y = 0, Ay ~ 1.875) as function of impact 
parameter b. 
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Figure 15: Comparison of our result for (^2) with STAR collaboration data[^. Open 
circles for RG and filled circles for QGP EoS. Whereas in the data points error bars are 
shown, in our estimates we have put 6v2. 
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b [fm] 


EoS 


event 


(^2) 


6V2 


3.5 


RG 


68 


0.0219 


0.0141 


QGP 


558 


0.0242 


0.0107 


7.0 


RG 


68 


0.0406 


0.0165 


QGP 


72 


0.0412 


0.0140 


10.0 


RG 


166 


0.0497 


0.0280 


QGP 


180 


0.0515 


0.0199 


12.0 


RG 


95 


0.0535 


0.0323 


QGP 


119 


0.0514 


0.0237 



Table 1: The average value and dispersions of the directed(i'i) and elliptic(i'2) flow coef- 
flcients for Au + Au coUisions at energy ^/s — 130 A GeV. 



b [fm] 


EoS 


event 


f 


6f 


3.5 


RG 


68 


0.229 


0.0024 


QGP 


55 


0.214 


0.0014 


7.0 


RG 


55 


0.231 


0.0032 


QGP 


58 


0.215 


0.0021 


10.0 


RG 


90 


0.233 


0.0041 


QGP 


119 


0.214 


0.0026 


12.0 


RG 


79 


0.234 


0.0047 


QGP 


100 


0.213 


0.0033 



Table 2: The slope parameter f and its dispersion 6f. 
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b [fm] 


EoS 


event 


A?/= 1.875 


A?/=3.00 


(n) 


dn 


dn/ (n) 


(n) 


dn 


6n/ (n) 


3.5 


RG 


68 


1026.5 


50.2 


0.049 


1619.2 


73.2 


0.045 


QGP 


55 


1557.5 


73.2 


0.047 


2548.4 


114.4 


0.045 


7.0 


RG 


55 


613.3 


49.5 


0.081 


977.4 


71.6 


0.073 


QGP 


58 


926.1 


81.1 


0.087 


1530.5 


123.7 


0.081 


10.0 


RG 


166 


312.8 


43.0 


0.137 


506.1 


65.8 


0.130 


QGP 


180 


437.5 


66.5 


0.151 


740.7 


103.9 


0.140 


12.0 


RG 


79 


162.8 


35.6 


0.219 


268.9 


56.2 


0.209 


QGP 


100 


220.1 


52.8 


0.240 


379.8 


85.2 


0.224 



Table 3: EoS dependence of the multiplicity fluctuation in two different Ay around y — 0, 
as function of the impact parameter b. 
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